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Abstract. - In this Letter we study the Kondo effect in an ultrasmall metallic grain, i.e. 
small enough to have a discrete energy-level spectrum, by calculating the susceptibility x of 
the magnetic impurity. Our quantum Monte Carlo simulations, and analytic solution of a 
simple model, show that the behavior changes dramatically depending on whether the number 
of electrons in the grain is even or odd. We suggest that the measurements of x provide an 
effective experimental way of probing the grain's number parity. 



Small metallic particles, typically of dimensions of few nanometers, show properties which 
are intermediate between atomic and bulk condensed matter systems. Because of the finite 
dimensions, the quasi-continuous spectrum in the bulk splits into discrete energy levels. When 
the average level spacing 5 is comparable with the other energy scales of the problem interest- 
ing new effects appear. For example the physical observables are parity dependent, i.e. differ 
if the number of electrons in the grain is even or odd. While metallic particles have been 
investigated in ensembles of grains in the past (see and references therein), only recently 
spectroscopy of single grains 2 has been achieved experimentally. In the case of supercon- 
ducting, ferromagnetic and/or antiferromagnetic systems, the effect of a finite level spacing 
is to enhance quantum fluctuations and consequently suppress the many-body condensation. 
As a result a superconductor will not possess a fully developed gap j3J, the magnetization of 
a ferromagnet will not be macroscopically large although distinct features reminiscent of 
the macroscopic order can be detected even in nanosize systems. 

Recently Thimm et al. [5], in their work on the so-called Kondo box, studied the effect 
of magnetic impurities in a small metal particle. They found that the Kondo resonance is 
strongly affected when the mean level spacing in the grain becomes larger than the Kondo 
temperature [SJ|B1[7]. The Kondo box shows the parity effect as well. The Kondo effect |S| 
is a paradigm in correlated electron systems. It originates from the interaction between 
the localized spin of a single magnetic impurity and the free electrons of the metal. In bulk 
systems Kondo correlations manifest in an increase of the resistance for decreasing temperature 
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T, above a characteristic scale Tk, the Kondo temperature. If, instead, one looks at the 
magnetic properties of metals hosting Kondo impurities, one effect is the rapid saturation, 
for decreasing T < Tx, of the impurity magnetic susceptibility x (measuring the effect of an 
external magnetic field on the impurity magnetization M) 0. Mesoscopic physics is another 
arena in which the Kondo effect is playing a major role. By now, several years since the 
theoretical prediction a series of beautiful experiments have seen it in different devices 

(see the reviews In these cases the resonant level is realized by suitable fabricated 

nanostructure and the conduction electrons are provided by the metallic electrodes. 

In this Letter we are interested in the Kondo box, and we want to address the parity effect 
by means of the analysis of the impurity susceptibility x(^)i since x(T) is very sensitive to the 
number of electrons in the system as demonstrated in superconducting dots |14ll5j . Detection 
of parity effects in various thermodynamical quantities, including the spin susceptibility, were 
measured for an ensemble of small, normal metallic grains |16| . In similar fashion parity effects 
can be measured in an ensemble of small grains containing Kondo impurities. As it will be 
discussed in more details, also in this case, due to parity effects, grains with odd number of 
electron will give the dominant contribution to x(T). 

The system is described by the Anderson model Hamiltonian 

H = ^{<^n ~ fJ,)cl ,^Cn,a + ^{f^d ~ fJ.)dlda + Und,^nd^ l + y ^(4,^d(T + Cn,adl) , (1) 

where cjj ^ {cn,cr) is the rising (lowering) operator in the grain electron state with spin a = ±1/2 
(or t, i) and energy e„ = ep + n6, labeled by the quantum number n starting from the lowest 
level Co < and dj. (d^) is the rising (lowering) operator in the impurity electron state with 
spin CT and energy ed with the associated number operator for spin up (down) in the impurity 
state Ud,^ = djdf {nd,i = rfjc^x)- The Fermi level is taken as the zero of the energy. The 
number of electrons N in the grain is selected by fixing eq, S, and the chemical potential to 
/X = for an odd number and fj, = S/2 for an even number, with (eg — /^) + [N — i)S/2 = 0. 
Only the levels up to — eo > are relevant to the problem and the density of grain electrons in 
the bandwidth 2|eo| is assumed constant. The energy U > accounts for the on-site repulsion 
in the localized state. The fourth term in Eq.QJ describes the hopping of an electron from a 
grain level e„ into the impurity level ed and vice versa, with a hybridization energy V. The 
level width T scales as F = ttV'^ /6 and remains finite for increasing S, because, as discussed in 
Ref. 0, both S and scales as the inverse of the grain volume. This is consistent with the 
increase of the probability (^ V"^) of finding an electron on the impurity site for decreasing 
grain volume. 

In order to put in context our analysis, let us recall some basic concepts concerning the 
physics associated with the Anderson model. The impurity site may be in one of the four 
different states: empty, singly occupied with up (down) spin, and doubly occupied. The 
magnetic properties of the impurity depend on the occupation probability of each state. As 
discussed in Ref. ^Jj, depending on T, one has different regimes. To illustrate how these 
different regimes develop, it is useful to consider first the case with no hybridization {V = 0) 
and ed < (Fig.^). For ksT ^ U+ed (where /cs is the Boltzmann constant), the four impurity 
states have the same probability to occur. Therefore the probability to have two electrons in 
the impurity level is P2 = 1/4, the average number of electrons in the impurity level is Nd = 1 
and the average of the squared value of the total spin in the impurity level is Af^ = 1/8. The 
impurity susceptibility as a function of T is, in this case, Xo{T) — Tk/{2T) (hereafter we use 
the dimensionless susceptibility xo = 4xfcB?V/(<?A*B)^j where g is the Lande factor and /is 
is the Bohr magneton), as in the case of a free orbital. For U + ed > ksT > —ed the double 
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Fig. 1 - Schematic representation of the energy levels for a system composed of a magnetic impurity 
embedded in a nanoscopic metallic grain in the Anderson model with V — 0. The possible states for 
the impurity level are: (a) empty level {Nd — 0) with M — 0; (b) or (c) singly-occupied level (A'^^ = 1) 
with M = — 1/2 (or spin I) or Af = 1/2 (or spin |); (d) doubly occupied level {Nd = 2) with electrons 
with anti-parallel spins (||) and M — 0. The corresponding energies are represented, in the bottom 
part, by bold lines for the states with the electrons filling the grain levels £„ up to the Fermi level 
and the impurity level empty (a), the impurity level singly-occupied with energy gain < (6 or 
c), the impurity level doubly-occupied by electrons repelling each other with energy U > (d). Thin 
lines represents the excited states with one or more grain electrons promoted to the first empty state 
above the Fermi level. 



occupancy is excluded, resulting in P2 ~ 0, TVd =i 2/3, AP ~ 1/6 and xo{T) ~ 2Tk/{'ST). For 
ksT < —Erf only the single occupied states have a finite probability to occur, with P2 ~ 0, 
Nd — 1, Af^ ~ 1/4 and Xo(r) ~ Tk/T. The impurity properties in this case are like those of 
a local (magnetic) moment^ i.e. of an isolated electron able to flip its spin a. 

If > 0, the electrons can tunnel between the impurity level and the grain levels. The 
virtual transitions to doubly occupied and empty impurity states lead to an effective antifer- 
romagnetic interaction J^f / between the grain electron spin density and the impurity spin. 

For a large grain's size {5 — s- 0), in first approximation (to the lowest-order of the per- 
turbation theory is Jeff ^ T/U. The higher orders can be taken into account by the 
renormalization group (RG) approach via an energy-dependent interaction Jg/ / (E) |51ll7j . 
The RG procedure shows that, at low T, Jeff{E) increases {strong coupling regime) ^Hl- Be- 
low Tfc, the impurity spin becomes effectively bound into a singlet (total spin zero) with the 
grain electrons' spin. As a consequence in the strong-coupling regime, the impurity suscepti- 
bility saturates, for T ^ 0, to xo{5 = 0, T = 0) ~ 0.4128 . 

In ultrasmall grains, a finite d introduces a further low-energy scale beside Tk- In order 
to describe the effect of finite level spacing on the Kondo resonance we use a quantum Monte 
Carlo approach proposed by Hirsh and Fye |19I2U| . According to this method, after integrating 
out the grain level degrees of freedom, the resulting impurity problem is mapped into a 
chain of two-states auxiliary variables (fictitious Ising spins), independent on each other, but 
interacting with an effective magnetic field, proportional to the impurity magnetization M. 
The length L of the chain is proportional to the time-evolution of the system. The method 
is exact in the limit of an infinite chain of Ising spins and gives results reliable within an 
approximation ~ (TL)"^. To study the properties of the impurity spin we calculate Xi P2, 
and as a function of T, and L (verifying, as a check, the exact relation 4M^ = Nd — 2P2 
valid at any T). We consider the symmetric case (i.e. with = —U/2). We estimate Tk for 
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5 — > 0, by fitting our calculations with the low-T approximation |S] 

4xfcB ^ 0-68 
{g^iB? T+IAU2TK 

where Tx is the only free parameter. In this limit our results (Fig. [21 recover, at fc^T ~ 
-ed = U/2, XoT/Tk ^ 1 (local moment regime) and, for T < Tk, Xo = 0.4128 (strong 
coupling regime). For S > 2kBTK, Xo (Fig. [2), P2 (Fig. OJ) and (Fig. ^ reveal a clear 
parity effect. The main features are the following. 

For an odd number N of electrons in the system (dotted lines in Fig. I2I4|I . i) the local 
moment regime extends down to the lowest T, with xo ^/T (Curie law); ii) xqT/Tk 
saturates to a value that increases with S, approaching the value 1 expected for F = 0, and iii) 
for T < T* and increasing S, P2 decreases toward as expected in the local moment regime. 

For an even N (solid lines in Fig. I2I4|I . iv) the Kondo effect appears to be enhanced and 
X saturates to a value that decreases with increasing 5. This can be seen as an increase of the 
effective Kondo temperature Tk{5) on the base of the Eq.|5J). v) For T < T* and increasing 
(5, P2 increases toward 1/4 as in the free orbital regime. 

Finally, vi) for ksT < V, with the hybridization energy V ^ 5^^^, the non-monotonic 
behavior of P2 is depressed when N is odd and emphasized when N is even; vii) for any iV, 
the onset of the parity effect increases with S and, empirically, the onset is at fc^T* ~ <5/5 for 
all the cases considered. The proportionality of the onset T* to S, is expected, because, for 
ksT > (5, the approximation of the grain levels with a band is valid. 

To understand the odd case, consider the limit S 00, equivalent to a system with an 
impurity level at energy ed < and a doubly occupied grain level at (zero) Fermi energy. The 
only possible process in this case is a transition of a grain electron to the impurity level, but 
for ksT < U this transition is suppressed, because the double occupancy in the impurity level 
costs an energy U. Therefore, the impurity electron behaves as an isolated magnetic moment, 
with P2 — and following the Curie law xo ~ ^/T (Fig. EJ, as in i). If S is finite, also the 
transition of the impurity electron above the Fermi level takes place, but requires an energy 
\ed\ + S and occurs with probability exp[—d/{kBT)]. Therefore, for large 5 and low T, this 
transition is suppressed and xqT/Tk — > 1 (inset Fig.jS)), as in the local moment regime, and 
P2 (Fig.lSJ, as in iii). For decreasing S, the system enters the local moment regime only 
a.t T < T* ^ 6, with xqT/Tk saturating at the value reached at T*. Therefore, the lower S, 
the lower the asymptotic value (inset Fig.l^J, as in ii). 

In the even case, when there is only one electron in the Fermi level, transitions between 
the Fermi level and the impurity level are allowed and this results in a strong hybridization 
between the two levels. We perform in this case an exact calculation of Xo{T) for S 00 
(Fig. 1^. In this limit, only the grain level at the Fermi energy couples with the impurity, 
resulting in an interacting two-level system. For very large, but finite S, this two-level system 
is still a good approximation since the grain levels away from the Fermi level can be taken into 
account in perturbation theory up to second order. Their main effect is to renormalize the 
impurity parameters. In particular, the energy of the impurity level is renormalized to a higher 
energy, closer to the Fermi level. This is just what is left of the logarithmic renormalization, 
familiar in the bulk system. The two-level system has six possible different many-body two- 
electron states: (1) with one electron with spin "f in the Fermi level and one with spin J, in the 
impurity |1) = c|(i||0); (2) with inverted spins with respect to the first state |2) = cj(i||0); (3) 
with two electrons with opposite spins in the Fermi level |3) = c|cj|0); (4) with two electrons 
with opposite spins in the impurity level |4) = d|(i||0); (5) with one electron with spin 1 in the 
Fermi level and one with spin t in the impurity |5) — ctd||0), and (6) with inverted spins with 
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Fig. 2 - The parity effect on the magnetic susceptibility. The (dimensionless) susceptibiUty xo, as 
function of the (dimensionless) temperature T/Tk, for a magnetic impurity embedded in a non- 
magnetic metalhc grain, shows a clear parity effect when the level spacing 5 is large: dotted lines 
correspond to odd number N of electrons in the system, solid lines to even A'^. We show the quantum 
Monte Carlo calculations for L = 300 for the symmetric Anderson model with T /{nksTK) — 1.62 <^ 
U/iksTK) ~ 20.45 and |eo| = 2556.25fcsTK, for S/iksTK) ~ 0.01 (bold line), recovering the 5 ^ 
limit, and, starting from the bold line and going outward, for 5/{kBTK) — 2.04 (rhombus) 3.37 (full 
circles), 5.11 (up triangles), 10.22 (open circles), 25.56 (crosses), and 255.57 (down triangles). The 
corresponding values of range from 2508 (solid line with rhombus) to 21 (dotted line with down 
triangles). Where not shown, errors are smaller than symbols. Inset: the same results times T /Tk 
are shown to emphasize the local moment regime for ksT ~ U /2. 




Fig. 3 Fig. 4 



Fig. 3 - The parity effect on the probability P2 of finding two electrons localized on the impurity. 
Symbols and lines are like those in Fig. |5| The calculations for 5 ~ O.OlfcsTx (bold line) show a 
non-monotonic behavior with a maximum for T < Tk- Dotted lines are for odd A'^; solid lines for 
even A''. The calculations are for L = 200. Comparison with calculations for L = 300 do not show 
any relevant finite-size effect for 5 > l.^ksTK- 

Fig. 4 - The parity effect on the squared impurity magnetization A/^. Symbols and lines are like 
those in Fig.|31 Since Nd^AM^ + 2P2, all our calculations have Nd — 1. 
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Fig. 5 - The exact susceptibility xo for even N and 5 —* oo, as a function of T /Tk- The calculations 
for the symmetric Anderson model with U — 20 AAksTK and V/ {ksTK) ~ 1.0 (solid line), 2.3 (dashed 
line) and 3.6 (long-dashed line), are shown in linear scale, in a) to emphasize the maxima of xo, and 
in double logarithmic scale, in b), to emphasize the high-T behavior xo ~ 1/r (the linear part in b). 



respect to the fifth state |6) — c|(i||0). These six states form three singlet (M = 0) and one 
triplet (Af = 0,±1) states. The singlet states are given by (|1) - |2))/\/2, |3), and |4). The 
triplet components are instead given by (|1) + |2))/a/2, |5), and |6). To probe the impurity 
susceptibility we couple exclusively the impurity spin with an external magnetic field h. The 
model, in this case, can be described by a diagonal block matrix, with a 4 x 4 block Ha, for 
the states from (1) to (4), and a 2 x 2 block Hb, for the states (5) and (6), with 
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2ed + C/ 



and H,= [ ^/^^ ) , (3) 



where the diagonal elements are the energies in each of the six states, and the square of an off- 
diagonal element is the transition probability between the corresponding states. By computing 
the energies Ei corresponding to the eigenstate of these two matrices, we calculate the impurity 
magnetization ■m{T,h) = exp[i?i/(fcBT)], where Mi is the impurity magnetization of 

the eigenstate i, and x(^) = idm/dh)h=o- 

We find that x ^ 1/^ for hgT > V (Fig.^). Indeed, for keT > V the states |2), 
1 5) and |6) have the same energy (for h = 0), while states |3) and |4) have higher energy, 
therefore the impurity level is singly occupied (local moment regime). For ksT < V, we 
find X saturating to a constant value, decreasing with increasing V (Fig. Ej). Indeed, as a 
consequence of the hybridization V, the ground state is given by the lowest-energy singlet, 
that has an energy difference with the triplet states larger than V. Therefore, x is constant 
and the asymptotic value is approximately given by the Curie law at ksT ~ V. Therefore, 
the larger V (or S), the smaller the saturation value of x for T ^ 0, as in iv). All these results 
are indeed consistent with our quantum Monte Carlo simulations (Fig. ^ , by keeping in mind 
that V scales as 6^^^ for fixed transition rate F. 
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Moreover, the actual diagonalization shows that the transition from the local moment 
to the saturation regime occurs via a non-monotonic behavior (with a maximum depending 
on the hybridization V). This result is not clearly detectable in our quantum Monte Carlo 
calculations for x (Fig- 13 j but is consistent with the non-monotonic behavior of P2 (Fig.|31l, 
more evident in the even case, that shows that the local moment regime extends to lower T, as 
in vi). In particular, to understand in the even case the behavior of P2 for finite decreasing S, 
we should consider more grain levels in our exact solution. For example, by including one grain 
level below the Fermi energy and one above, the many-body states are four-electron states, 
due to three grain levels and one impurity level. As noted before, the hybridization favors a 
singlet ground state. In this case the weight of the singlet state with doubly-occupied impurity 
level is decreased, because the singlet state manifold has a larger number of components. As 
a consequence, the probability P2 decreases with decreasing 6 (Fig. |3J), as in v). 

In summary, we have shown that, a magnetic impurity embedded in a small metallic grain, 
has a striking different behavior depending on the parity of the electron number in the grain. 
In particular, the Kondo effect is strongly enhanced if the number of electrons is even, while 
it is strongly depressed if the number is odd. Therefore, the measurement of x of a magnetic 
impurity embedded in an ultrasmall metallic grain is an effective way of probing the number 
parity of the electrons in the grain. 
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